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ABSTRACT 


This  paper  reports  sensitivity  improvement  of  a 
stationary  GPS  receiver  using  software  GPS 
approach.  The  experiments  were  conducted 
using  both  simulated  and  real  GPS  signals.  The 
sensitivity  improvement  can  be  as  much  as  1 0 
dB.  This  paper  discusses  the  software  algorithms 
developed  to  perform  the  code  acquisition  and 
tracking  that  enabled  the  sensitivity 
improvement.  The  commercial  GPS  receiver  has 
an  unusual  sampling  frequency  of  5.455657x1 06 
Hz.  The  paper  will  discuss  techniques  used  to 
handle  such  odd  sampling  rate  in  the 
softwaealgorithms,  also,  the  determination  of  the 
time  skew(fme  time)  between  clock  phase  and 
initial  phase  of  C/A  code  when  working  with  this 
odd  sampling  rate.  The  techniques  discussed  in 
this  paper  is  not  only  limited  to  work  with  the 
specific  commercial  receiver.  They  can  be 
applied  to  software  receiver  with  any  sampling 
frequency.  Experimental  set  up  and  procedures 
used  to  evaluate  the  sensitivity  improvements  are 
also  discussed. 

Introduction 

There  are  several  GPS  manufactures  that  provide 
GPS  OEM(Original  Equipment  Manufacture) 
modules  with  a  RF  front  end,  A/D  converters, 
channels  of  correlator,  and  a  postprocessor. 
These  modules  can  be  connected  to  a  portable  or 
a  desk  top  computer.  The  computer  can  process 
information  provided  by  the  GPS  module  for  a 
variety  of  applications.  Due  to  the  ever 
increasing  computational  power  of  personal 
computers,  for  those  applications  that  require  the 
use  of  a  personal  computer,  we  may  take  one 
step  further  in  exploiting  the  use  of  the 
computer.  That  is,  we  can  simplify  the  GPS 
module  and  move  GPS  receiver  functions  to  the 
PC  to  implement  a  software  GPS  receiver.  This 
can  be  achieved  by  using  the  RF  front  end  chips 
on  these  OEM  modules  to  down  convert  the 
input  signal  to  the  IF  (Intermediate  Frequency) 
band  and  then  digitize  the  signal.  The  rest  of  a 
GPS  receiver  functions  can  all  be  accomplished 
by  software  on  the  computer.  Not  only  does  this 
approach  require  less  hardware  components,  it 
also  has  better  performance  in  the  code  initial 
phase  acquisition  because  the  frequency  domain 
block  processing  can  be  used.  The  tracking  loop 
can  also  be  more  adaptive  to  signal  strength, 
noise  condition,  and  application  dynamics. 


Figure  1.  Hardware  setup  for  receiver  sensitivity 
evaluation. 

SENSITIVITY  COMPARISON 

The  front  end  used  in  this  study  is  a  part  of 
Motorola  Model  Ml 2  Oncore  GPS  receiver.  The 
sampling  frequency  is  5.455657xl06  Hz.  We 
modified  this  front-end  so  that  IF  output  could  be 
tapped  and  digitized.  The  experimental  set-up  is 
shown  in  Figure  1.  The  top  portion  of  Figure  1 
shows  the  setup  of  experimenting  with  simulator 
data.  A  hardware  GPS  simulator  (Model 
STR4790  from  Global  Simulation  System,  Inc.) 
was  connected  to  the  input  of  the  receiver.  Eight 
satellite  signals  were  programmed,  and  the  signal 
strength  in  five  of  them  are  very  close.  The 
power  level  of  the  signal  could  be  varied  by  an 
attenuator.  We  found  that  the  Motorola  receiver 
stopped  tracking  signals  when  the  attenuator  was 
set  at  71  dB.  The  receiver  could  not  be  cold 
started  until  the  attenuator  setting  was  reduced  to 
61  dB.  Digitized  signals  were  collected  with  the 
attenuator  settings  from  61  to  73  dB.  The 
software  approach  could  acquire  and  track  the 
signal  at  72-73  dB,  demonstrating  that  the 
sensitivity  of  the  receiver  was  improved  by  about 
1 0  dB  while  the  same  hardware  front  end  was 
used. 

The  bottom  portion  of  Figure  1  shows  the  setup 
of  working  with  real  data.  A  special  front-end 
arrangement  was  made  for  this  purpose.  A 
passive  antenna  followed  by  an  attenuator  was 
placed  before  the  first  RF  amplifier.  With  this 
set-up,  when  the  attenuator  was  changed  one  dB, 
the  signal-to-noise  ratio  of  the  receiver  is  also 
changed  one  dB.  Therefore,  the  amount  of 
improvement  could  be  determined.  Similar 
results  were  obtained  as  in  the  case  of  using  the 
simulator.  The  minimum  attenuator  setting  was 
recorded  before  the  Motorola  receiver  could  cool 
start.  From  this  minimum  setting,  the  attenuation 
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(1) 


was  increased  to  1 0  dB  while  the  software 
approach  could  still  acquire  and  track  the  signals. 
Thus,  the  sensitivity  improvement  of  software 
receiver  was  about  1 0  dB  when  working  with 
real  data. 

C/A  CODE  INITIAL  PHASE  AND 
FREQUENCY  ACQUISITION 

In  the  software  approach,  an  initial  processing  is 
to  further  down  convert  the  digitized  data  from 
IF to  base  band.  The  base  band  frequency  can 
vary  ±5KHz  due  to  satellite  Doppler  shifts  even 
when  GPS  receiver  is  stationary.  It  is  well 
understood  that  one  ms  of  coherent  data 
processing  can  provide  1  kHz  frequency 
resolution  while  10  ms  of  coherent  data 
processing  can  provide  1 00  Hz  frequency 
resolution.  This  paper  uses  one  ms  coherent  data 
processing  to  perform  the  circular  correlation. 
Ten  consecutive  1-ms  coherent  correlation 
results  from  the  same  time  bin  can  be  used  to 
perform  10-point  FFT(Fast  Fourier  Transform) 
to  coherently  integrate  10  ms  of  data[1]. 


l-ms  C/A  code  in  signal 
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Figure  2.  Signal  C/A  code  and  data  block 
alignment 


Frequency  in  KHz 


Figure  3  C/A  code  spectrum  (with  carrier) 

Since  the  sampling  frequency  is  5. 45 5657MHz, 
the  number  of  samples  in  one  period  of  Doppler 
free  C/A  code  is  5455.657.  If  Doppler  shift  is 
present,  the  number  of  samples  in  one  period  of 
C/A  code  will  change  by  this  amount: 


f 

m  =  5455. 657( — — — ) 

f  l  f  d 

where  fi  =  1575.42  MHz  is  the  LI  frequency  and 
fd  is  the  Doppler  frequency.  According  to  (1), 
Doppler  frequency  of  ±5  KHz  correspond  to  m 
values  of  5455.657±0.0173  point.  It  is  not  a 
significant  change.  Therefore,  the  number  of  data 
in  one  C/A  code  period  is  very  close  to 
5455.657.  But  it  is  not  an  integer.  In  order  to 
perform  the  block  circular  correlation,  we  have 
to  make  data  block  size  an  integer.  As  a  result, 
5456  data  points  are  taken  as  a  block  for 
frequency  domain  processing.  For  every 
millisecond,  5456  data  are  read  in  and  down- 
converted.  As  shown  in  Figure  2,  assuming  the 
first  point  of  the  first  data  block  is  coincidently 
lined  up  with  the  initial  phase  of  the  C/A  code  in 
input  signal.  The  first  point  in  the  second  data 
block  is  0.343  point  off  the  initial  phase  of  the 
second  C/A  code  in  the  signal.  If  we  continue 
taking  data  in  this  way,  the  first  point  of  the  third 
data  block  will  be  0.686  point  off  the  initial 
phase  of  the  third  C/A  code  in  the  signal.  This 
offset  is  greater  than  0.5  point.  To  ensure  the 
offset  to  be  less  than  half  of  a  sample  interval, 
the  first  point  of  the  third  data  block  is  taken 
from  the  last  point  of  the  previous  block.  Now, 
the  first  point  of  the  third  data  block  is  -0.314 
point  off  the  initial  phase  of  third  C/A  code  in 
the  signal.  The  following  data  blocks  are  taken 
based  on  the  same  principle  to  ensure  the  first 
point  of  a  data  block  and  the  initial  phase  of  C/A 
code  of  the  signal  are  within  ±0.5  sample  point. 
After  data  are  taken,  a  5456-point  FFT  (Fast 
Fourier  Transform)  is  applied  to  them.  Figure  3 
is  a  FFT  spectrum  of  C/A  code  with  IF  at 
1 .364Mhz.  It  shows  that  less  than  half  of  the 
spectrum  contains  significant  information. 
Therefore,  only  half  of  the  spectrum  is  kept.  The 
other  half  is  casted  out  since  it  contains  half  of 
the  total  noise  and  only  a  small  portion  of  the 
useful  informational  Keeping  only  half  of  the 
spectrum  also  reduces  the  computation  time.  The 
local  C/A  code  reference  is  also  digitized  into 
5456  points  and  applied  the  FFT.  The  same  half 
of  the  complex  conjugate  spectrum  is  stored  for 
repetitive  use.  The  acquisition  uses  10  ms 
coherent  integration  and  then  uses  as  many  as 
non-coherent  integrations  needed  to  acquire  the 
initial  phase  of  C/A  code  and  the  Doppler 
frequency. 

Due  to  the  Doppler  effect,  the  location  of  the 
initial  phase  of  C/A  code  in  each  data  block  is 
different.  In  order  to  non-coherently  integrate 
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tens  of  blocks  of  data  together,  the  C/A  code  of 
the  input  data  in  each  block  has  to  line  up  with 
the  initial  phase  of  C/A  code  in  the  first  block. 
Time  shift  needed  for  each  block  of  data  can  be 
computed  as  following: 

rD  =  Y  SD  =  vup(-jlrfn  X  Td)  (2) 

TD  is  the  time  shift  needed  and  fD  is  the  Doppler 
frequency.  T  is  time  delay  between  the  initial 
time  of  the  current  data  block  and  the  initial  time 
of  the  first  block  of  integration.  SD  is  the 
frequency  domain  equivalent  of  TD .  /n  is  the 
frequency  of  kernel  function  of  the  FFT.  For 
easier  Matlab  implementation,  this  time  shift  is 
handled  in  the  frequency  domain.  To  obtain 
correlation  result  of  a  1  ms  data  block,  a  circular 
correlation  is  performed.  This  is  done  by 
performing  the  inverse  FFT  of  the  product  of  SD, 
the  spectrum  of  the  data  block,  and  the  stored 
complex  conjugate  of  the  C/A  code  spectrum. 
After  10  sets  of  circular  correlation  are 
performed,  the  results  of  the  same  time  bin  are 
coherently  integrated  through  a  10-point  FFT. 
There  are  2728  10-point  FFT  results  after  this 
process.  They  form  a  2728x10  time  frequency 
spectrum  matrix.  This  matrix  covers  the  Doppler 
frequency  range  of  ±500  Hz.  The  detailed 
approaches  have  been  discussed  in  a  previous 
paper[2].  In  order  to  cover  ±500+nxl000  Hz,  the 
input  data  have  to  be  multiplied  by  exp(-j27i 
nlOOOt).  The  result  is  then  applied  to  FFT  to  get 
new  spectrum  and  undergo  the  same  process 
mentioned  above.  Another  approach  is  to  rotate 
the  spectrum  of  input  data  left  for  n  bins.  In  order 
to  cover  all  possible  Doppler  frequencies,  as 
many  as  1 1  (n=  -5,-4, ..,0,.. ,4, 5)  matrices  have  to 
be  created.  The  result  is  a  2728x10x11  matrix. 
However,  if  Doppler  of  satellite  can  be 
estimated  through  the  ephemeris  data  from 
external  sources,  the  size  of  the  matrix  can  be 
greatly  reduced. 


Figure  4.  Filter  shape  for  Fourier  Transform 
based  filter  bands  (or  bins) 


Figure  5.  Definition  of  the  correspondent  bin 


Figure  6.  The  sensitivity  of  the  original  matrix. 


Figure  7.  The  sensitivity  after  the  deficiency  fix 
between  bands 


Frequency  in  H  z 


Figure  8  The  sensitivity  after  all  deficiency  fix 

There  is  one  deficiency  associated  with  the 
above  approach  and  needs  compensation 
techniques  to  overcome  it.  Lets  use  filer  bands  as 
a  series  of  1  KHz  filters  and  use  filter  bins  as  a 
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Incoherent  integration  gain 


series  of  100  Hz  filters.  Figure  4  shows  the  filter 
shape  for  the  Fourier  transformation  based  filter 
bands  (or  bins).  If  a  signal  with  the  Doppler 
frequency  falls  between  the  frequency  bins,  the 
signal  power  will  be  split  into  two  bins  and 
phases  of  the  components  in  two  bins  are  almost 
inverted  when  integration  length  as  large  as  1 0 
data  points  or  more.  The  splitting  of  the  signal 
into  two  separate  bins  leads  to  3.92  dB 
sensitivity  degradation.  Figure  4  can  also 
represent  the  filter  shape  between  two  1  KHz 
bands.  Therefore,  the  receiver  will  experience 
3.92  dB  loss  if  this  Doppler  frequency  is 
between  two  1kHz  bands.  As  a  result,  if  a 
signal’s  Doppler  frequency  happens  to  fall 
between  frequency  bins  and  between  two  1kHz 
bands,  the  loss  can  be  7.84  dB.  The  combined 
result  is  shown  in  Figure  6. 

Our  software  receiver  uses  an  efficient  way  to 
reduce  the  loss  between  the  1kHz  bands.  This  is 
achieved  by  creating  a  2728x10x10  time- 
frequency  spectrum  matrix.  The  matrix  is 
obtained  by  computing  the  difference  between 
the  same  two  bins  in  filter  bank  m  and  m+1, 
scaled  by  a  factor  of  1.414.  Figure  5  can  be  used 
to  explain  how  this  matrix  is  obtained.  The  first 
5  bins  of  bank  1  and  last  5  bins  of  bank  1 1  have 
no  correspondent  bins.  They  are  ignored  from 
now  on.  The  scale  factor  is  used  to  counter  the 
increase  of  the  noise  power  (3  dB)  during  the 
differencing  operation.  If  the  signal  frequency 
falls  between  bands,  the  improvement  of 
sensitivity  is  3  dB  using  this  new  time-frequency 
spectrum  matrix.  But  the  sensitivity  for  signals 
with  the  frequencies  at  the  center  of  the  band  will 
be  reduced  by  3  dB.  This  is  if  a  signal  frequency 
falls  at  the  center  of  band  m,  then  the  signal  level 
at  band  m+1  is  near  0,  resulting  in  no 
improvement  in  the  differencing  signal  power. 
The  noise  power,  however,  will  be  increased  by 
3  dB  during  the  differencing  operation. 

Therefore,  the  signal  to  noise  ration  will  be 
degraded  by  3dB.  In  order  to  avoid  this  3  dB 
loss,  the  differencing  and  scaling  operation  is 
only  applied  to  the  5  bins  that  falls  between  two 
adjacent  bands  of  the  original  spectrum  matrix. 
Figure  7  shows  the  result  obtained  from  the 
matrix  created  in  this  way. 


Figure  9.  Gain  improvement  verse  Non-coherent 
integration  number 

In  our  implementation,  No  second  matrix  is 
created.  We  only  use  the  above  principle  to 
modify  all  5  bins  falling  between  bands  in  the 
original  matrix. 

In  order  to  reduce  the  loss  between  bins,  a  new 
modified  2728x10x10  time-frequency  spectrum 
matrix  is  created  by  subtracting  each  bin  with  its’ 
adjacent  bin  and  scaling  with  1.414.  Now  we 
have  two  2728x10x10  time  frequency  spectrum 
matrix  waiting  for  following  10-ms  coherent 
integration  result  to  non-coherently  add  on.  After 
all  required  integrations,  the  maximum  of  two 
matrices  are  detected  and  its’  time  and  frequency 
are  selected  as  initial  phase  of  C/A  code  and 
course  Doppler  frequency  which  has  accuracy 
within  ±50  Hz.  Figure  8  shows  the  result  of  the 
sensitivity  improvement  combining  two  matrix. 
The  number  of  non-coherent  integration  is 
determined  based  on  the  desired  sensitivity  of 
the  receiver.  Figure  9  is  used  to  determine  the 
number[6]. 

Because  we  want  to  use  one  second  coherent 
integration  for  weak  signal  tracking,  the  Doppler 
frequency  has  to  be  accurate  within  ±0.5  Hz. 
Since  the  initial  phase  of  the  C/A  code  and 
coarse  carrier  frequency  have  been  determined, 
the  carrier  frequency  of  data  in  each  block  can  be 
down  converted  to  ±50  Hz.  The  C/A  code  of  the 
signal  can  be  stripped  off  by  multiplying  the 
down-converted  signal  with  local  reference  C/A 
code.  One  block  of  data  can  be  integrated  into 
one  point.  After  one  second,  1 000  integrated  data 
can  be  obtained.  By  performing  FFT  on  the 
squared  integrated  data,  we  can  achieve  Doppler 
frequency  accuracy  of  within  ±0.5  Hz.  To  begin 
weak  signal  tracking,  we  also  need  to  determine 
the  initial  phase  of  navigation  code.  The 
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algorithm  to  find  the  initial  phase  of  navigation 
code  has  been  discussed  in  a  previous  paper[3]. 

LOW  DYNAMIC  AND  HIGH  SENSITIVITY 
C/A  CODE  AND  FREQUENCY  TRACKING 


initial  phase  offset  is  defined  as  the  offset 
between  the  initial  phase  of  C/A  code  and  the 
first  sample  of  data  in  the  data  block.  The  initial 
phase  offset  measurement  will  be  discussed  later. 

*  r"  Offset  x 


Input  data  sampleLS^^ 


|  f0f  0;'  | 


|  .] 


Tim  er 
pointer 


1-ms  Doppler  free  C/A 
code 


5456  |  545(/ 


1-ms  data  block  for  a  specific 
satellite 


5456  | 

5456 


Initial  phase  of  signal  C/A  code 


First  sample  of  data  block 


First  sample  of  reference  data  block 


Initial  phase  of  reference  C/A  code 


Figure  10.  Satellite  specific  data  block  for 
tracking 

In  order  to  maintain  the  phase  continuity  of  the 
carrier,  5456  points  of  data  have  to  be  read  in 
continuously.  There  is  a  pointer  and  a  mismatch 
accumulator  associated  with  each  satellite  signal. 
Assuming  there  is  no  Doppler  shift,  and  the 
initial  phase  of  the  C/A  code  of  a  satellite  in  the 
signal  is  1000  points  off  from  the  beginning  of 
the  first  data  block.  The  pointer  is  pointing  at 
1000.  This  is  shown  in  Figure  10.  The  mismatch 
between  1  period  of  code  length  in  the  signal  and 
the  data  block  is  0.343  sample.  After  two  blocks, 
the  mismatch  is  0.686  which  is  greater  than  0.5 
point.  The  pointer  is  updated  to  999  which  is  the 
closest  integer  point  to  the  real  initial  phase  of 
code  in  the  signal  and  the  mismatch  accumulator 
is  updated  to  -0.314.  Let’s  also  look  at  a 
different  scenario.  The  initial  phase  of  the  C/A 
code  of  a  satellite  in  the  signal  is  1  point  off  the 
beginning  of  the  first  data  block  and  the  pointer 
update  is  needed.  The  pointer  is  supposed  to  be 
updated  to  0,  but  there  is  no  0  point  in  the  block. 
Under  this  circumstance,  the  pointer  is  updated 
to  5456  and  the  pointer  is  tracking  the  initial 
phase  of  C/A  code  of  the  following  period  .  In 
the  program,  the  pointer  and  the  mismatch 
accumulator  are  constantly  updated  to  maintain 
the  mismatch  between  the  timer  pointer  and  the 
real  initial  phase  of  code  in  the  signal  to  within 
±0.5  point.  The  mismatch  due  to  Doppler  is 
important  in  this  process  and  can  not  be  ignored 
because  the  mismatch  update  is  an  accumulative 
process.  Every  second  a  5456-point  digitized 
Doppler  free  C/A  code  with  a  proper  initial 
phase  offset  is  generated  as  local  reference  to 
‘perfectly’  strip  off  C/A  code  in  the  signal  in  the 
following  second.  As  shown  in  Figure  1 1,  the 


Figure  11.  Definition  of  offset 

The  high  sensitivity  tracking  algorithm  uses  1- 
second  coherent  integration.  The  integration 
starts  at  the  initial  phase  of  any  navigation  bit. 
Each  millisecond,  the  data  is  read  in  and 
multiplied  by  a  locally  generated  C/A  code  and  a 
locally  generated  complex  IF  carrier  signal 
(within  1  Hz  accuracy).  The  multiplication 
results  are  added  together  to  generate  one 
complex  data  point.  Twenty  1-ms  data  are  added 
together  to  generate  one  navigation  data.  The 
algorithm  to  decode  the  weak  navigation  data  has 
been  discuss  in  a  previous  paper[4].  After  the 
navigation  data  is  decoded,  the  navigation  data 
are  multiplied  with  the  navigation  code  and  then 
added  to  produce  1 -second  coherent  integration 
result.  The  amplitude  of  the  1 -second  coherent 
integration  is  called  the  1 -second  prompt  gate 
result.  The  same  procedure  can  be  used  to 
produce  early  and  late  gate  results. 

The  early  gate  result  is  generated  by  using  C/A 
code  in  the  signal  that  is  one  sample  ahead  of  the 
locally  generated  C/A  code  reference,  while  the 
late  gate  result  uses  the  signal  that  is  one  sample 
later. 
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Figure  12  Late  gate  satellite  data  and  C/A 
reference  alignment. 
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Figure  13  Late  gate  results  from  figure  12  in 
correlation  function 


For  the  ease  of  illustration,  the  local  reference  is 
generated  without  offset  in  Figure  12.  Figure  13 
shows  the  late  gate  results  of  the  first  three 
milliseconds  on  the  correlation  triangle  function. 
They  are  located  around  1+x  and  based  on  the 
mismatches,  d’s,  which  we  have  discussed 
before.  The  averaged  late  gate  correlation  result 
is  (Yl+Y2+Y3)/3  at  x+7+  (~o.314+0+0.343)/3. 
From  this  figure,  we  can  derive  the  1000-point 
averaged  late  gate  result  and  its’  location  as 


1000 


IX  at 


1  +  X  + 
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IX  (3) 

n= 1 


if  we  denote 

X  —  x  H - V  d  (a\ 

1000  ( j 

The  point  (Xa+\  ,  Yal)  is  still  on  the  right  side  of 
the  triangle.  By  the  same  approach,  we  can 
derive  that  the  early  gate  result  will  be  (Xa-l  , 

Yai )  and  also  on  the  right  side  of  the  triangle. 
Therefore,  Xa  can  be  calculated  by  using 
following  method[5]. 


Yea  Ypa  Yu 


Xa-d  Xa  Xa  +d 


Figure  14.  Correlation  function  and  early, 
prompt,  and  late  gate 


r  „  i  -  x  .  -  d 

Y  ae  1  +  X  a  ~  d 

( 1  -  r  X  1  -  d  ) 


(5) 


Here  Yae  and  Yai  are  the  1-sec  early  and  late  gates 
respectively  and  d  and  Xa  are  in  unit  of  chip.  For 
our  case,  d  is  equal  to  183.3/977.5.  After 
calculation,  Xa  has  to  be  converted  back  to  unit 
of  sample.  The  C/A  code  initial  phase  offset  x 
can  be  derived  from 


2  1000 

x  =  X - V  d  (in  unit  of  sample)  (6) 

1000 

This  offset  is  important  for  the  pseudo  range 
calculation  and  the  generation  of  C/A  code 
reference  for  the  next  second  of  data.  This 
completes  the  C/A  code  tracking  process. 


In  frequency  domain,  a  50-point  FFT  are  applied 
to  the  navigation-code-stripped  complex 
navigation  data  as  mentioned  above.  By 
comparing  the  amplitude  of  the  peak  and  the 
amplitude  of  its  two  neighbor  frequency  bins,  the 
frequency  offset  can  be  calculated[4].  The 
frequency  of  the  peak  and  frequency  offset  are 
added  to  the  IF  carrier  frequency  to  generate  the 
locally  generated  IF  carrier  signal  for  next 
second.  This  completes  the  frequency  tracking 
process. 


TRACKING  PERFORMANCE 
EVALUATION 


8  Seconds  of  Simulated  data  were  used  to 
evaluate  the  tracking  performance.  At  each  C/N0, 
1 00  runs  with  different  noise  distributions  are 
performed.  The  results  are  listed  in  Table  1.  The 
software  receiver  can  track  the  C/A  code  and 
frequency  but  navigation  code  starts  appear  error 
after  C/N0  <24dB.  Since  the  number  of  runs  in  is 
rather  small,  it  is  difficult  to  draw  a  firm 
conclusion. 
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Table  1  Tracking  performance 


C/N0  (dB) 

Mean 

error(ns) 

Standard 

deviation(ns) 

25 

-1.1 

14 

24 

-2.1 

16 

23 

-2.9 

18 

22 

-7.5 

20 

21 

-6.5 

22 

20 

-3.2 

27 

19 

1.3 

30 

18 

4.3 

35 

CONCLUSION 

This  software  receiver  has  demonstrated  a  1  OdB 
sensitivity  improvement  for  stationary  receivers. 
This  means  that  in  a  very  weak  signal 
environment,  it  can  acquire  satellite  C/A  code 
and  carrier  frequency  and  smoothly  transfer  from 
acquisition  mode  into  tracking  mode.  This 
receiver  is  not  limited  to  weak  signal 
environment.  With  very  minor  software  code 
changes,  a  phase  lock  loop  can  be  implemented 
for  strong  signals  under  high  dynamic 
conditions.  This  receiver  with  the  carrier  phase 
tracking  loop  can  track  2  minutes  of  a  hardware 
simulator  generated  data  with  GPS  receiver  on 
an  airplane  pulling  8  G  S-turn  at  the  speed  of 
400m/s.  The  details  of  the  high-dynamic 
application  are  beyond  the  scope  of  this  paper. 

REFERENCE 

(1)  DAVID  M  LIN  and  JAMES  B.Y.  TSUI 
“Direct  P(Y)-Code  Acquisition  Algorithm 
for  Software  GPS  Receivers”,  Proceedings 
of  ION  GPS  99”  September  14-17,  1999, 
Part  1,  pp.  363-368 

(2)  DAVID  M  LIN  and  JAMES  B.Y.  TSUI 
“Acquisition  Schemes  for  Software  GPS 
Receiver”,  Proceedings  of  ION  GPS  98” 
September  15-18,  1998,  Part  1,  pp.  317-326 

(3)  DAVID  M  LIN  and  JAMES  B.  Y.  TSUI  “ 

A  Software  Receiver  for  Weak  Signals”, 
IEEE  International  Microwave  Symposium, 
May  20-25,  2001, 


(4)  DAVID  M  LIN  and  JAMES  B.Y.  TSUI  “A 
Weak  Signal  Tracking  Technique  for  a 
Stand-Alone  Software  GPS  Receiver”, 
Proceedings  of  ION  GPS  02”  September  24- 
27,  2002 

(5)  JAMES  B.  Y.  TSUI  “  Fundamentals  of 
Global  Positioning  System  Receivers”,  John 
Wiley  &  Sons,  INC,  pp  182-183 

(6)  Lin,  D.  M.,  Tsui,  J.  B.  Y.,  Liou,  L.,  and 
Morton,  Y.  T.,  “Sensitivity  limit  of  a  stand¬ 
alone  GPS  receiver”,  Proceedings  of  2002 
ION  GPS  conference,  Portland,  Oregon, 
Sept.  2002. 


8 


